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1. Introduction 

The goal of this paper is to investigate a fairly general enumeration problem related to 
the theory of knots, links and tangles: we want to count objects which live in 3-dimensional 
space and are (loosely) made of a certain collection of "ropes" , some of which open (with 
fixed endpoints) and some closed on themselves, intertwined together in an alternating 
way. As usual in knot theory, these objects will be considered up to topological equivalence 
(deformation or ambient isotopy), and represented by their projections on the plane; we 
shall then classify them according to the (minimal) number of crossings, the number of 
connected components, and the way the various "external legs" connect to each other, see 
for example Fig. 1. 




Fig. 1: A diagram with 3 open and 1 closed line, intertwined together with 
7 crossings. 

Without going into too much detail for now, we see that a convenient way to keep 
track of the number of connected components and of the connections of the external legs 
is to use colors, see Fig. 2. 




Fig. 2: Coloring the diagram of Fig. 1. Open lines have fixed colors (distinct 
from each other), whereas closed lines have arbitrary color. 
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The colors allow us to distinguish the various external legs and add an extra power 
series variable in the theory (the number of colors n) to count separately objects with 
different numbers of connected components. 

The idea to use colors was already suggested in [1,2] and present in the work [3]. 
At this stage it is natural to define a matrix model whose Feynman diagram expansion 
will produce such diagrams with n colors. Here we shall give a unified quantum field- 
theoretic treatment of this 0(n)-invariant matrix model, which simplifies and generalizes 
the equations obtained in [3] (section 2 below). In particular, it gives a practical way to 
do the enumeration by computer; this procedure was recently used in the numerical work 
[4]. 

Even though the matrix model we propose is fairly natural, since as we shall see it 
is the most general quartic 0(n)-invariant matrix model with a single trace in the action, 
it is in general unsolvable (or at least unsolved). It can be thought of as describing a 
statistical model on random dynamical lattices; more precisely, it is a model of fully packed 
loops drawn in n colors on random tetravalent planar diagrams with weights attached to 
vertices (intersections or tangencies of loops). Even the corresponding model on a regular 
(flat) square lattice is not fully understood. However it is tempting to speculate on its 
universality class; and that putting it on random lattices will correspond to the usual 
coupling of two-dimensional conformal field theory to gravity, which allows to predict 
the critical exponents of the theory based on the KPZ relation [5]. This in turn leads to 
various conjectures on the asymptotic number of large links and tangles, made in [2], which 
have been checked numerically in [4]. We shall not come back to these conjectures here, 
but instead produce exact analytic solutions of two particular cases of our matrix model 
(section 3 below): the classical case n = 1 (no colors), with some generalizations of the 
results of [1]; and the case n = —2, which is interesting because its asymptotic behavior 
cannot be obviously guessed by the universality arguments mentioned above. 

2. General principle 

We assume the reader familiar with the concept of links and tangles. Let us recall here 
that once projected on a plane, they give rise to planar diagrams with tetravalent vertices 
which must be "decorated" to distinguish under/over-crossings. Link diagrams are closed, 
whereas tangle diagrams have external legs. The diagrams are said to be alternating if 
one meets undercrossings and overcrossings alternatingly as one follows the various closed 
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loops of the diagram. The alternating property allows to ignore the decorations of the 
vertices since they can be recovered from the diagram alone (up to a mirror symmetry for 
the closed diagrams, see below). 



2.1. Definition of the 0{n) matrix model 

As in [2] and [3], we start with the following matrix integral over N x N hermitean 
matrices 

Z<«)(n,,) ^ ff[ dM„e^'' ^ i ^^»^^'*^"*^') (2.1) 

a=l 

where n is (for now) a positive integer. The integral is normalized so that Z^^\n^ 0) = 1. 
The partition function (2.1) displays a 0{n) symmetry where the Ma form a vector of 
0{n). 

Expanding in power series in g generates Feynman diagrams with double edges ( "fat 
graphs") drawn in n colors in such a way that colors cross each other at the vertices. By 
taking the large limit one selects the planar diagrams/ which are closely related to 
alternating link diagrams, cf Fig. 3. 




Fig. 3: A planar Feynman diagram of (2.1) and the corresponding alternat- 
ing link diagram. 



More precisely, the large A^ "free energy" 



logZ(^)(n,ff) 

F{n,g)= hva — 2.2 



^ Note that the next orders in the 1/N expansion of the free energy log Z''^^ would correspond 
to link diagrams drawn on thickened surfaces of higher genus, cf [6]. 
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is a double generating function of the number /^.p of alternating link diagrams with k 
connected components and n crossings (weighted by the inverse of their symmetry factor, 
and with mirror images identified): 

oo oo 

^(^'^?) = EE^;p^V (2.3) 

k=lp=l 

Note that it is clearly possible to analytically continue Z^^^ (n, g) to arbitrary values of 
n (using, for example, a Hubbard-Stratonovitch transformation) so that Eq. (2.3) still 
holds. In particular, the counting of knot diagrams is given by Fi-p and can be obtained 
by formally taking the limit n ^ 0, in the spirit of the replica method. Also, if n is an even 
negative integer one can write fermionic analogues of (2.1), see section 3.2, which display 
Sp{\n\) symmetry. 

If one is interested in counting objects with a weight of 1, one cannot consider the 
free energy which corresponds to closed diagrams, but instead correlation functions of the 
model which generate diagrams with external legs: these are essentially tangle diagrams. 
Typically, we shall be interested in the two-point function 

where the measure on the the matrices Ma is given by Eq. (2.1) and a is any fixed index, 
which generates tangle diagrams with two external legs; and the connected four-point 
functions 

ri(n,^)= hm /i-tr(M,M6)A (2.5.1) 

JV^oo \ iV / 

r2{n,g) = lun^ (^^ir{MlMl)J - G{n,gf (2.5.2) 

where a and h are two distinct indices, which generate tangle diagrams with four external 
legs of type 1 and 2 (see Fig. 4). Note that the freedom to replace a link diagram with 
its mirror image by inverting all under/over-crossings is, in the case of correlation func- 
tions, removed by fixing conventionally the first crossing encountered starting from a given 
external leg. 

Let us briefiy mention for now that the definition of G{n,g) again assumes n to be 
a positive integer, and has a natural continuation to any n; however the definitions of 
ri{n,g) are only meaningful for n integer greater or equal to 2, and there is a difficulty 
associated to this, which will be explained in section 2.3. 
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Fig. 4: Tangles of types 1 and 2. 
2.2. Renormalization of the 0{n) model 

The model presented above is not sufficient to properly count colored tangles. Es- 
sentially, this comes from the fact that there is not a one-to-one correspondence between 
diagrams and the objects they are obtained from by projection. This generates a redun- 
dancy in the counting since to a given knot will correspond many (an infinity of) diagrams, 
each counted once. In the case of alternating diagrams one can distinguish two steps to re- 
move this redundancy. First one must find a way to restrict ourselves to reduced diagrams 
which contain no irrelevant crossings (Fig. 5a)); such diagrams will have minimum num- 
ber of crossings. It turns out to be convenient to introduce at this point a closely related 
notion: a link is said to be prime if it cannot be decomposed into two pieces in the way 
depicted on Fig. 5 b). It is clear that at the level of diagrams, forbidding decompositions 
of the type of Fig. 5 b) automatically implies that the diagram is reduced; and we shall 
therefore restrict ourselves to prime links and tangles. 




Fig. 5: a) An irrelevant crossing, b) A no n- prime link. 

There may still be several reduced diagrams corresponding to the same link: according 
to the flyping conjecture, proved in [7] , two such diagrams are related by a finite sequence 
of flypes, see Fig. 6. 




Fig. 6: A fiype. 
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To summarize, there are two problems: a) the diagrams generated by applying Feyn- 
man rules are not necessarily reduced or prime; b) several reduced diagrams may corre- 
spond to the same knot due to the flyping equivalence. A study of Figs. 5 and 6 shows 
that this "overcounting" is local in the diagrams in the sense that problem a) is related 
to the existence of sub-diagrams with 2 external legs, whereas problem b) is related to a 
certain class of sub-diagrams with 4 external legs. Clearly such graphs can be cancelled 
by the inclusion of appropriate counterterms in the action. We are therefore led to the 
conclusion that we must renormalize the quadratic and quartic interactions of (2.1). Now 
renormalization theory tell us that we should include in the action from the start every 
term compatible with the symmetries of the model, since they will be generated dynam- 
ically by the renormalization. In order to preserve connectedness we only look for terms 
of the form of a single trace. A key observation is that, while there is only one such 
quadratic 0(n) -invariant term, there are two quartic 0(n)-invariant terms, which leads to 
a generalized model with 3 coupling constants in the action (bare coupling constants): 

'{n,t,gi,g2)= [[dMae V 

a=l 

(2.6) 

The Feynman rules of this model now allow loops of different colors to "avoid" each other, 
which one can imagine as tangencies (Fig. 7). 



Fig. 7: Vertices of the generalized 0{n) matrix model. 

We define again the correlation functions G{n,t, gi, §2) and ri{n,t, gi, §2) (Eqs. (2.4) 
and (2.5)), and want to extract from them the counting of colored alternating tangles with 
external legs. 

The idea is to find the expressions of t{g), gi{g) and g2{g) as a function of the renor- 
malized coupling constant g, in such a way that the overcounting is suppressed and the 
correlation functions are generating series in g of the number of colored tangles. At leading 
order, we shall have t{g) = 1 -|- o(l), gi{g) = g + o{g) and (72 ((7) = o{g) so that we recover 
the original model (2.1). However there will be higher order corrections which correspond 
to the counterterms. 
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Let us consider ^(^f) first. It is clear that one must remove all two-legged subdiagrams, 
that is impose 

G'(n,t,(7i,^72) = l (2.7) 




+ 



+ 



Fig. 8: Decomposition of the two-point function. Reexpanding in powers of 
t — 1 will cancel the powers of S iff t = 1 -|- E. 

Let us see more explicitly how this fixes t{g). Noting that (Fig. 8) 

G 1 



t 



i2.i 



where E is the generating function of IPI (one-particle irreducible, i.e. which cannot be 
made disconnected by removing one edge) two-legged diagrams, one finds equivalently that 



tig) = 1 + Eig) 



(2.9) 



i.e. the counterterms generated by t{g) must cancel all IPI two-legged subdiagrams. This 
is almost a tautology; notice however that one must not cancel all two-legged subdiagrams, 
since one-particle reducible diagrams would be subtracted multiple times. 






Fig. 9: Breaking a flype into two elementary flypes. 

Next, we must consider the flyping equivalence. Again, it is important to notice that 
a flype can be made of several "elementary" flypes (Fig. 9), an elementary flype being by 
definition one that cannot be decomposed any more in this way. In the terminology of 
QFT, an elementary flype consists precisely of one simple vertex connected by two edges 
to a non-trivial H-2PI (two-particle irreducible in the horizontal channel) tangle diagram. 
Non-trivial means not reduced to a single vertex; H-2PI means that the tangle diagram 
cannot be cut into two pieces containing the left and right external legs respectively, by 
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removing two edges. We therefore need to introduce auxiliary generating functions H[{g), 
H2{g) and l^'lfi*) non-trivial H-2PI tangles of type 1, of type 2 and of type 2 rotated 
by 7r/2 respectively. Only these must be included in the counterterms. It is now a simple 
matter to consider all possible insertions of elementary flypes as tangle sub-diagrams of a 
diagram; taking into account the two types of tangle sub-diagrams and the two channels 
(horizontal and vertical), we find (Fig. 10) that the renormalization of gi and g2 is simply: 

gi{g) = g{l-2H!,ig)) (2.10.1) 
g2i9) = -9{H[ig) + VM) (2.10.2) 







8i K^bK^b 8 K^bK^t 




S2 KKM^M, 





Fig. 10: Counterterms needed to cancel flypes. 



All that is left is to find the expressions of the auxiliary generating functions in terms 
of known quantities. They are easily obtained by decomposing the four-point functions in 
the horizontal and vertical channels, and will not be rederived here (the reader is referred 
to e.g. [3] for details). 
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Ho ih H\ =1 — 

" ' (iT9)(i + A±ri 
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(i-(7)(i-F(n + i)r2 + ri) 
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(2.11a) 
(2.11b) 



2.3. Summary and discussion 

We shall now summarize and rewrite more explicitly the formulae found previously, 
as well as discuss their implications. 

Let us assume that for a certain n, we have computed the free energy F{n,t, gi, g2). 
What can we extract from the formulae of the previous paragraph, and how? 

First, let us differentiate F; we find 

d F 

G = -2^- (2.12) 
ot n 

as well as two other quantities. 



F =4^- = - hm (Itr V(M„M6)M (2.13.1) 
n N—>oo \ iV ^ — ' / 

\ a,b I 



d F 1 

'1 = 

ogi n 



d F 1 

Fo = 2 



- lim ( ^tr Vm^MH (2.13.2) 

\ a,b I 



dg2 n 

According to the equations of motion, these three quantities are not independent: 

tG = 1 + giFi + 2g2F2 (2.14) 

Comparing (2.13) with the definition (2.5) of the Tj, we see that there are two different 
choices of basis of the four-point functions;^ using 0(n)-invariance of the measure it is easy 
to relate them: 

Fi = uFi + 2{r2 + C^) (2.15.1) 
F2 = A + (n + l)(r2 + G^) (2.15.2) 

These relations also have a simple diagrammatic interpretation, which proves in particular 
that they are valid for any (complex) n. One observes that relations (2.15) can be inverted 
to extract Fi and A only if n 7^ 1, —2. These two cases will be the object of study of 
the next section, they are the first in the series of bosonic / fermionic matrix models; and 
as will be shown these are the values of n for which the model possesses only one quartic 
0(n)-invariant, contrary to the generic case. For now we simply observe that for n = 1, 



^ Using the Fi as the preferred basis is not only natural diagramatically; it is also imposed by 
the structure of the equations such as (2.10) and (2.11). 
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Fi = F2 and therefore according to (2.13), the free energy F is a function of gi +2g2 only; 
while for n = — 2, Fi = —2Fi and F is a function of gi — (72 only. 

Once we have computed G, A and 12, we can slightly simplify the renormalization 
equations using the obvious scaling property: 

G(n, t, gi, g2) = ^G{n, 1, gi/t^, ^72 A') (2.16a) 

ri{n,t,g^,g2) = ^r,(n, 1, V*', W*') (2-166) 
Combining this with Eq. (2.7) results in fixing t{g): 

t{g) = Gin, l,gi{g)/t{g)\g2{g)/t{gf) (2.17) 

At this stage the three unknowns t{g), gi{g), g2{g) only appear through the combinations 
hi{g) = gi{g) / t{gY , ^2(5') = 92{g)/t{g)'^; in particular we have the following expressions 
for the Fi = Fi{n,t{g), gi{g), g2{g)): 

_ r,in,lJn{g)J,2ig)) 
' G{n,l,h,{g),h2{g)r 

We only need to solve the two remaining renormalization equations (2.10), which we rewrite 
here: 

h,{g) Gin, 1, hi{g),h2{g))^ = g{l - 2H^ig)) (2.19.1) 
h2{g)G{n,l,hi{g),h2{g)f = -g{H[{g) + VM) (2-19.2) 

where the auxiliary generating functions are still given in terms of the Fi by Eqs. (2.11). 

Finally, solving Eqs. (2.19) gives access to the Fi, which are the generating series of the 
numbers of prime alternating tangles of type i. However, we can go further. By computing 
other correlation functions in the model and composing them with the solutions t{g), gi{g), 
g2{g) of the equations above, one can extract the generating functions of the number of 
alternating tangles with an arbitrary number of external legs. The correlation functions 
we consider are traces of non-commutative words in the Ma of degree 2k (for 2k external 
legs). We usually restrict ourselves to connected correlation functions (free cumulants in 
the language of free probabilities) , which exclude configurations in which some strings have 
no crossings with the other strings and can be pulled out altogether. This choice is only a 
matter of taste. 
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(2.18) 



For example, there are five 0(n)-invariants of degree 6, except, as before, for a finite 
set of values of n for which there are fewer: only 4 for n = —4, 3 for n = 2, 2 for n = —2, 
1 for n = 1. They are given by: 



:i = lim ( —tv{MaMbMcMaMbM^) 



l2 = lim ( —MMaMbM^MaM^Mb) 



h = lim ( —tr{MaMaMbMcMbMc 



14 = lim ( —MMaMMMaM^M, 



15 = lim ( —tr{MaMaMbMtM,M, 



— disc, terms 

— disc, terms 

— disc, terms 

— disc, terms 

— disc, terms 



(2.20.1) 
(2.20.2) 
(2.20.3) 
(2.20.4) 
(2.20.5) 



(a, 6, c distinct) and give rise to the various six-legged diagrams depicted on Fig. 11. 








Fig. 11: The five types of tangles with 6 external legs. 
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3. Application: two solvable cases 

There are currently two values of n for which the corresponding matrix model has 
been exactly solved: n = 1 and n = 2. The case n = 1 is particularly important since 
it corresponds to counting all alternating tangles regardless of the number of connected 
components; we shall investigate it here in detail, generalizing known results [8,1]. 

The application of the 0{n = 2) matrix model (also known as six- vertex model on 
dynamical random lattices) to knot theory has already been made in [3], using slightly 
different methods than in the present paper, and we shall not come back to it. 

However, we have found earlier that aside from n = 1, there is another special value 
of n, namely —2, for which a simplification in the model occurs and we can expect some 
exact analytic results. We shall present below an analysis of this case. 

3.1. The case n — 1: the usual tangles, and more 

As an illustration of the general principle developed above, we present an elementary 
solution of the case n = 1, that is the counting of alternating tangles. Since there are 
no colors one cannot distinguish the way the various external legs are connected; the 
correlation functions available to us will be specified by the number of external legs only. 

Note that this solution, which generalizes the original calculation of the number of 
prime alternating tangles with 4 external legs found in [8], is technically different from it. 

We start by setting n = 1 in the definition of the partition function (Eq. (2.6)); we 



where qq = gi + 2g2. The fact that the partition function only depends on a particular 
combination of gi and g2 is consistent with what was found in Section 2.3 and related to 
the existence of only one quartic 0(n) -invariant for n — 1. The most general "planar" 
correlation functions of the model are of the form 



find: 




(3.1) 




(3.2) 



for which we introduce the generating function: 




(3.3) 
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and the corresponding connected correlation functions G^i-, whose generating function is 
the inverse function A(a;): 

= - (3-4) 

Among them we have the two-point function G = G2 = G2 and the connected four-point 
function F = G% — G4 — 2G2 which is nothing but the generating function of all tangles 
(regardless of type): F = F1+2F2. Since we do not have access to A and F2 separately, we 
need to recombine the equations of Section 2.2 so that only F appears in them. Fortunately, 
this turns out to be possible; taking (2.10.1)-|-2x (2.10.2) results in 

go{9)^9{l-2H{g)) (3.5) 

where H{g) = H'2{g) + H[{g) -\- V2{g) is the generating function of all H-2PI non-trivial 
tangles. Eq. (3.5) can of course be derived directly in a manner similar to Eqs. (2.10), by 
simply disregarding the types of the tangles i.e. how the outgoing strings are connected to 
each other inside the tangle. 

Setting n = 1 in Eq. (2.116), we also find that 

so that for n = 1 (and n = 1 only) we have a closed subset of equations. 

We now turn to the solution of our matrix model. We do not repeat the calculation 
of the G21 here since it is a standard result of matrix models, see [9]. Starting from the 
following expression: 



(A) = \{t\- goX^ -{t- goX" - goA)VX''-2A^ (3.7) 



with A = ^!zV^E}^ solution of 
3 go 



3A^go - 2At + 4 = (3.8) 



we find that 



In particular G = |^(4 — ^), and since G = 1 according to Eq. (2.7), we can express t as 
a function of A: 

t=^i2A-3) (3.10) 
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Similarly, using the explicit expression of F = G4 — 2, plugging it into Eqs. (3.5), 
(3.6), and using Eqs. (3.8), (3.10) to express t and go in terms of A results in the following 
fifth degree equation for A: 

32{l-g)-64{l-g)A + 32{l-g)A^-4{l + 2g-g^)A^ + Qg{l-g)A^-g{l-g)A^ = (3.11) 

A{g), specified by Eq. (3.11) and A{g = 0) = 2, is a well-defined analytic function of g in 
a neighborhood of ^ = 0. The data of A{g) is enough to recover all correlation functions 
since we have, combining Eqs. (3.9) and (3.10): 

G2i = 2A^-^^^^^^{3e-{e-l)A) (3.12) 



(£ + 2) 



Similarly, one can extract the connected correlation functions, using the fact that A(u;) 
satisfies a cubic equation (cf Eq. (3.7)); after a tedious calculation one finds 

G^, = |(A - 2)^-1 {3i-2-{i- 1)A) (3.13) 

where q is a constant (which already appeared in [9]): 

This concludes the calculation of the generating series of the number of tangles with any 
given number of external legs. In the appendix, the first few orders of G4, Gq, Gg are 
given. 

Let us now discuss the asymptotic behavior of the coefficients of the various series 
for which we found an exact expression. All are simple polynomials in A{g), so that we 
need to study the latter only. As can be easily checked, the singularity of A{g) closest to 
the origin is the usual singularity of 2D pure gravity, which is given by goc = 4/27 and 
tc = 4/3, so that Ac = 3, and, plugging into Eq. (3.11), 

V2IQQI-IOI 
9c = (3.15) 

We expand A around g ] gc and find 



A = 3-aigc- g)'/^ + bigc-g) + 0{{gc - gf'^) (3.16) 

14 



with (a > 0) 



2 2877137 + 7087^21001 



= (3.17a) 

339696 ^ ' 

5(99397733 + 2127733a/2T00T) 

~ 901510722 

This provides the leading singular part of G^f 

Gli = reg + j^a[b + ^(£ - 2)a') {g, - gf/^ + • • • (3.18) 
which finally yields the large order behavior of C^f if (^2^ = Z^^i lU\p9^ then 



p— >CX) 



3 Q /, ,1 



[i - 2)! 



a(6+ ^(£-2)a2)p-5/2^3/2-f (3.19) 



For £ = 2 this result coincides with the theorem 1 of [8] . Note that for any i the asymptotic 
behavior is the same up to a constant. One can of course send i and p to infinity in a 
correlated manner to obtain a non-trivial scaling limit (here, i oc p^^^); but the result 
is known to be universal and is therefore the usual scaling loop function of pure gravity, 
which will not be reproduced here. 

3.2. The case n — —2: a fermionic matrix model 

For n negative even integer, it is natural, in the spirit of supersymmetry, to look for 
realizations of our model of colored links under the form of a fermionic matrix model with 
Sp{—n) symmetry. Let us show how this works in the simplest case, that is n = —2. 

Our fields will be a "complex fermionic matrix", that is a matrix ^ = (\I/ij) where 
the are independent Grassmann variables, and its formal adjoint = (^'ji), which 
together form the fundamental representation of Sp{2). We apply to them the usual 
rules of Berezin integration. We must next look for 5'p(2)-invariant quadratic and quartic 
invariants of the form trP(\E', \1/^). Since the \E' are non-commutative objects, one must 
consider arbitrary tensor products of the (dual of the) fundamental representation; however 
the trace property combined with the anti-commutativity of the matrix elements implies 
that an elementary circular permutation must have eigenvalue —1 in this representation. 
Very explicitly, there are one quadratic invariant, _^t^^ two independant quartic 
invariants, say ^r^rt^^t _ ^t^^^t _ $$t$t$ _|_ $t$$t$ and ^r^r^rt^t _ ^t^^^t _ 
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$$t$t$_l_$t$t$$. However it is clear that the first quartic invariant is stable by circular 
permuation and therefore its trace is zero. We are thus left with the following expression: 

It is no surprise that the partition function only depends on one coupling constant go, since 
the analysis of section 2.3 has shown us that for n = —2 all large quantities depend only 
on the combination (71 — (72; and indeed, by direct inspection one can identify go — gi — 92- 
As in the case n = 1, we have access to only one four-point fonction 

r = lim ( — tr**^*^* ) = A - (3.21) 

Again, a "miracle" happens in that a particular subset of the renormalization equations 
becomes closed for n = —2; namely, taking (2.10.1) — (2.10.2) 

^o(^) = 9{l - 2H'^{g) + H[{g) + Vi{g)) (3.22) 

a combination of the H-2PI diagrams appears, which can be related to F via the use of 
Eqs. (2.11): 

We now briefly describe how to compute integral (3.20) in the large N limit. We can 
set t = 1 without loss of generality, as explained in section 2.3. We perform the standard 
Hubbard-Stratonovitch transformation by introducing a hermitean matrix A: 

The gaussian integral over ^' and ^'^ can then be performed: 

Z(^) {t, go) = J det(l 1 + ^/=^(^ ® 1 + 1 ® ^)) e~f ^^^^ (3.25) 

We recognize at this point the usual 0{—2) fully packed (non-intersecting) loops model^ 
[10]. We perform the change of variables: M = {A — oq)^ with ao = — 1 which 
absorbs the determinant, resulting in 

Z(^)(t,^o) = y"dMe-f MVM + «o)' (3.26) 



^ Of course this might have been expected from the start, since for any n the 0(n) model of 
fully packed non-intersecting loops corresponds to the particular case g"! = of our model, and 
therefore here qq = —52- 
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One can then diagonalize M and compute the integral over eigenvalues using standard large 
N saddle point techniques. The resolvent of M is given by a complete elliptic integral of 
the third kind; in particular, 

where K and E are complete elliptic integrals of the first and second kind with modulus 
k, and the coupling constant go is given by 

^0 = -^^((2 - k'')K - 2E) (3.28) 

(cf also [11] for a similar solution). Finally, inserting the expression of = ^^"^2 + 
1 obtained from (3.27) into renormalization Eqs. (3.22)-(3.23) results in a (/-dependent 
transcendental equation for the modulus k!^ . This equation is too complicated to be solved 
exactly; however it can be easily solved numerically order by order. Finally, F is the 
desired generating function Fi — F2 of tangles; in the appendix we present the first few 
orders of the expansion of F[g). 

We now turn to the asymptotic behavior of the coefficients of F{g). It is known 
that the 0{—2) matrix model does not have any critical point of the usual form of 2D 
quantum gravity; for example in [11], the solution of the 0{—2) model of dense loops, 
equivalent to ours, is studied in detail when the elliptic modulus k is in the range [0, 1], 
which corresponds to ^ro £ [—^i 0] in our notations, and no singularity is found. This does 
not mean, of course, that F{g) has no singularity, since ^ro (and g) can move in the whole 
complex plane. Generically, these singularities are square root type singularities and given 
by the equation ^ = 0. It turns out that this equation has plenty of solutions, though 
one can unfortunately not write them down analytically. Numerically, one finds that the 
the solutions with smallest modulus of g are: 

£fc ~ -0.239 ±0.135z (3.29) 

They are pairs of complex conjugate solutions: this indicates oscillatory behavior of the 
series, which is due to the fact that n < 0. We therefore find that ii F = Yl'^=i 1p9^^ 

-ip ~ Re(cst p~^^^g~P) (3.30) 

where est ^ —0.237 ± 0.090z. It would be interesting to find a physical interpretation of 
this critical point a la Yang-Lee. 
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Finally, let us note that one could combine the results of n = +2 [3] and of n = —2: 
this would give rise to a model of oriented tangles in which one counts separately tangles 
with odd and even numbers of connected components. Since the coefficients of r{g) (and 
presumably also of ri{g) and r2{g) separately) in the case n = —2 satisfy, according to 
Eq. (3.30), 7p = (9(3.64 . . .^) whereas in the case n = +2 these coefficients are of the order 
6.28 . . .P, we conclude that there are, up to exponentially small corrections, as many odd 
tangles as there are even tangles. . . 

Appendix A. Tables for n = 1 and n = —2 up to 32 crossings. 
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Tab. 1: Table of the number of prime alternating tangles (n = 1) with 4, 6, 
8 external legs. 
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Tab. 2: Table of the coefficients of T (n = —2 tangles with 4 external legs). 
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